Astronomy & Astrophysics manuscript no. 19410-newj-astroph 


©ESO 2013 


March 8, 2013 





Fragmentation in the Massive Star-Forming Region 

IRAS 19410+2336* 



J. A. Rodon 1 , H. Beuther 2 , and P. Schilke 3 



(N 

o 

C 
3 



< 

O 

Ok 
6 

CO 



m 
m 

O 

(N 



X 



1 European Southern Observatory, Alonso de Cordova 3107, Vitacura, Casilla 19001, Santiago 19, Chile, 
e-mail: jrodon@eso.org 

2 Max-Planck-Institut fur Astronomie, Konigstuhl 17, 69117 Heidelberg, Germany. 

I. Physikalisches Institut, Universitat zu Koln, Zulpicherstr. 77, 50937 Koln, Germany. 



Received ; accepted 



ABSTRACT 



Context. The Core Mass Functions of low-mass star-forming regions are found to resemble the shape of the Initial Mass Function 
(IMF). A similar result is observed for the dust clumps in high-mass star forming regions, although at spatial scales of clusters that 
do not resolve the substructure that is found in these massive clumps. The region IRAS 19410+2336 is one exception, having been 
observed at spatial scales on the order of ~ 2500 AU, sufficient to resolve the clump substructure into individual cores. 
Aims. We investigate the protostellar content of IRAS 19410+2336 at high spatial resolution at 1.4 mm, determining the temperature 
structure of the region and deriving its Core Mass Function. 

Methods. The massive star-forming region IRAS 19410+2336 was mapped with the PdBI (BCD configurations) at 1.4 mm and 3 mm 
in the continuum and several transitions of formaldehyde (H 2 CO) and methyl cyanide (CH 3 CN). The H 2 CO transitions were also 
observed with the IRAM 30 m Telescope. 

Results. We detected 26 continuum sources at 1 .4 mm with a spatial resolution down to ~ 2200 AU, several of them with counterparts 
at NIR and MIR wavelengths, distributed in two (proto)clusters. With the PdBI CH 3 CN and PdBI/IRAM 30 m H 2 CO emission we 
derived the temperature structure of the region, ranging from 35 to 90 K. Using these temperatures we calculated the core masses of 
the detected sources, ranging from ~ 0.7 to ~ 8 M Q . These masses were strongly affected by the spatial filtering of the interferometer, 
filtering out a common envelope with ~ 90% of the single-dish flux. Considering only the detected dense cores, and accounting for 
binning effects as well as cumulative distributions, we derived a Core Mass Function, with a power-law index /? = -2.3 ± 0.2. We 
resolve the Jeans length of the (proto)clusters by one order of magnitude, and only find little velocity dispersion between the different 
subsources. 

Conclusions. Since we cannot unambiguously differentiate protostellar and prestellar cores, the derived CMF is not prestellar. 
Furthermore, because of the large fraction of missing flux, we cannot establish a firm link between the CMF and the IMF. This 
implies that future high-mass CMF studies will require to complement the interferometer continuum data with the short spacing in- 
formation, a task suitable for ALMA. We note that the method of extracting temperatures using H 2 CO lines becomes less applicable 
when reaching the dense core scales of the interferometric observations because most of the H 2 CO appears to originate in the envelope 
structure. 
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1. Introduction 

Our understanding of the structure of the cold, dense interstel- 
lar medium (ISM) in star-forming regions has improved in the 
last years. In those regions the ISM exhibits a clumpy, often 
filamentary structure with density maxima at the sites of star 
formation. To represent this structure quantitatively we use the 
Core Mass Function (CMF). In this paper we will refer to "core" 
as the small (diameter D ~ 0.01 pc), dense condensation that 
will form individual stars or small multiple systems, while with 
"clump" we denote structures that may form (proto)clusters and 
may therefore be more massive and larger than cores. In this 
sense, cores may be considered as a subset of clumps. 

Sub-mm observati ons of low-mass star- forming regions 
such as Serpens (e.g., iTesti & Sargentl [T998I) . Orion B (e.g., 



* Based on observations carried out with the IRAM Plateau de Bure 
Interferometer and the IRAM 30m Telescope. IRAM is supported by 
INSU/CNRS (France), MPG (Germany) and IGN (Spain) 



iMotte et alj|2001|), Aqu ila (e.g jKonvves et alj|2010h and p Oph 
(e.g. JMotte et al J 19981). as well as near-infrared extinction maps 
(e.g., lAlves et all 120071) . show that their CMFs resemble the 
shape and intri nsic mass scale of the stella r initial mass func- 
tion (IMF; e.g., Salceter 1955tlKroupal2002l) . This suggests that 
these dense cores would be the immediate precursors of stars, 
and that by applying a more or less constant core-to-star mass 
conversion efficiency we can obtain the IMF from the CMF. 

In the case of Massive Star- Forming (MSF ) regio ns we 
have, for example, the analysis of iReid & Wilson! (120061) . They 
gathered the publ ished mass es of the MSF regio ns M8 (e.g, 
iTothill et al.ll2002h. M 17 (e g. jReid & Wilsonl2006l) NGC7538 
(e.g.. IReid & Wilson! 120051). W43 (e.g ., IMotte et al.ll2003h and 
RCW 106 (e.g., Mookerjea et alj |2004). tracing spatial scales of 
clumps that correspond to (proto)clusters rather than to individ- 
ual cores, and tested the fit of several functional forms for their 
clump mass functions. They found that in those cases, the best 
fit was obtained by a double power law, having a mean power- 
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law exponent for the high-mass end consistent with the Salpeter 
IMF. This would again imply that by an almost one-to-one 
mass conversion efficiency the IMF could be obtained from the 
clump mass function, as in the case for low-mass star-forming 
region s . Similar analysis w ere conducted by e.g., IShirlev et alj 
d2003l):lBeltran et al.ld2006l). F urthermore, a theoretical study by 
Chab rier & Hennebelle ( 20101) on the relationship between the 
CMF and the IMF suggests a tight correlation between the two. 
Once more, this implies that the IMF would be defined by the 
CMF in the early stages of evolution. 

However, this one-to-one relationship may not hold as, for 
example, some clumps must fragm ent to produce the obs erved 
quantity of multiple stellar systems dGoodwin et all 120071) . The 
relatively large distances (> 2 kpc) of most of the known MSF 
regions require a spatial resolution of about 1" to resolve the 
clumps into cores with sizes below ~ 0.1 pc. That resolution in 
the (sub)mm regime is only achievable with the interferometric 
technique. So far, only a few MSF regions have been observed 
in the (sub)mm with spatial resolutions good enough to resolve 

r ' II — ~i I II rn 

individual c ores (e . g.JBontemps et al.H2 010, Fontani et al.H2Q09L 
iRodon et al]|2008l iRathborne et al.ll2008l iBeuther et al.ll2006T) . 
and only for one source, IRAS 194 10+2336, has it been possi ble 
to determine a core mass function (IBeuther & SchilkeU2004l) . 

The young MS F region I RAS 1 9410+2336 is at a dis- 
tance of 2. 16 kpc dXu et all 120091) and has an integrated 
bolometric luminosity of about 1O 4 L . It is a very active 
star-forming site, with sources detected from X-rays down 
to ra dio wavelength s . It h a s H2O and CH3OH m aser emis- 
sion dSridharan et all 120021 IBeuther et all l2002dl) and X-ray 
2002X orma- 



sources (Beuther et al 



tion of intermediate-to-high mass stars. The region is embedded 
within a cluster of over 800 components detected at NIR wave- 
lengths dMartm-Hernandez et all 120081: lOiu et al.L 120081) . and it 
has a rich and energ etic outflow component with multiple out- 
flows detected in CO dBeuther et all I2002d.l2003h . IBeuther et alJ 
d2002bl) found that the large-scale mm emission shows two mas- 
sive gas clumps roughly aligned in a north-south direction that 
splits into several subsou rces with increasing spatial resolution 
dBeuther & Schilkdl2004l ). 

Wit h their studies of th e mm continuum at high spatial res- 
olution, IBeuther & Schilkd (12004 were able to derive the mass 
function of IRAS 19410+2336, resulting in a Salpeter-like dis- 
tribution. However, the strongest caveat in the derivation of that 
mass function was the fact that a uniform dust temperature was 
used in the calculation of the masses. Although they argue that 
the dust temperature distribution should not vary strongly, they 
also warn that changes in the temperature of the cores would 
result in a somewhat flattened slope. 

We have revisited IRAS 19410+2336 observing the mm con- 
tinuum at high-spatial resolution and obtaining molecular-line 
emission of known temperature tracers, to determine via several 
methods a temperature structure for it and in the end derive a 
more robust mass function. 



2. Observations 

2.1. Interferometric 

We observed the two protostellar clusters of IRAS 19410+2336 
with the PdBI in the B (Feb-2005), C (Dec-2004/Mar-2005) 
and D (Apr-2005) configurations, comprising baselines from 
20m to 330m. This translates into projected baselines rang- 
ing from ~ IkX to ~ 120 &/1 at 3 mm and from ~ 15&/1 to 
~ 240 kA at 1 .4 mm. The phase centers were set at RA( J2000) = 



19 /, 43 ffl 10.7 I ; Dec(J2000) = 23°44'58.4" for the "Northern" 
(proto)cluster and at RA(J2000) = 19 A 43 ffl l 1.2 s ; Dec(J2000) = 
23°44'03.2" for the "Southern" (proto)cluster. The continuum 
was mapped at 3 mm and 1 .4 mm. 

Table 1. Observed molecular transitions and rms of the respec- 
tive maps. 
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Spectral 
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'-'up 


rms 
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Transition 




(GHz) 


Resol. 


(K) 


(mjybeam ') 








(MHz) 




North 


South 


CH 3 CN(6 


-5 ) 


110.383 


0.18 


18.5 


10 


17 


CH 3 CN(6i 


-5,) 


110.381 


0.18 


25.7 


10 


17 


CH 3 CN(6 2 


-5 2 ) 


110.375 


0.18 


47.1 


10 


14 


CH 3 CN(6 3 


-5 3 ) 


110.364 


0.18 


82.9 


10 


14 


H 2 CO(3 , 3 


- 2o, 2 ) * 


218.222 


0.36 


21.0 


29 


36 


H 2 CO(3 2 ,2 


- 2 2 ,i) * 


218.476 


0.36 


68.1 


19 


30 


H 2 CO(3 2 ,i 


- 2 2l0 ) 


218.760 


0.36 


68.1 


23 


34 



Notes. {a> For a spectral resolution of 0.5 km s 1 . " J) rms values obtained 
after combining 30 m and PdBI data. 

The 3 mm receiver was tuned in the upper sideband and the 
1 .4 mm receiver in the lower sideband. With this spectral setup 
we observed the H2CO and CH3CN transitions described in 
TableQ]with a maximum spectral resolution of 0.5 km s" 1 , adopt- 
ing a systemic veloc i ty Vi sr = 22.4kms~ 1 dRidge & MooreL 
l200Tt rneftrunk et all [19981 

The phase and amplitude calibrators were 1923+210 and 
2023+336 and the flux calibrators were 3C273, 2200+420, 
1749+096 and 1741-038, adopting for the calibration the flux 
values from the SMA flux monitoring of these quasars^. The data 
were calibrated with CLIC and then imaged with MAPPING, 
both part of the GILDAS package^ The spectra were processed 
with CLASS90, also from the GILDAS package. 

After imaging and deconvolution, the resulting synthesized 
beams for the continuum data are 1.2" x 0.8" at 1.4 mm and 
2.2" x 1 .6" at 3 mm. At the given distance of 2.2 kpc that means 
a spatial resolution of ~ 2200 and ~ 4200 AU, respectively. The 
continuum data for the Northern (proto)cluster have rms noise 
levels cr ~ 0.8 mJy beairT 1 and cr ~ 0.4mJy bearrT 1 at 1.4mm 
and 3 mm respectively, while in the southern (proto)cluster 
the rms noise levels are cr ~ l.OmJy beairT 1 and cr ~ 
0.4mJybeam _I at 1.4mm and 3 mm respectively. For the line 
data, the synthesized beams and noise levels are detailed in 
Table □ 



2.2. Short spacings 

We obtained IRAM 30 m observations of H 2 CO (3 ,3 - 2 ,2) and 
H 2 CO(32,2 -2 2 ,i) towards IRAS 19410+2336 in November 
2007, taken in the on-the-fly mode with both HERA hetero- 
dyne receivers tuned at 218.222GHz. For the backend we used 
the VESPA correlator, assigning 2 of its spectral bands to each 
HERA receiver, with a channel spacing of 320 kHz and a band- 
width of 80 MHz, resulting in a spectral resolution of 0.5 km/s 
at 218 GHz. 

The data were processed with CLASS90. The single-dish 
uv-data were combined with the interferometric uv-data using 



http://smal.sma.hawaii.edu/callist/callist. html 
|http : //www . iram. fr/IRAMFR/GILDAS | 
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MAPPING. After imaging and deconvolution we obtained a syn- 
thesized beam of 1.6" x 1.0" for the combined data. The rms 
levels of the combined data are shown in Table Q] 

3. Results 

3.1. Millimetric Continuum 

We established a detection threshold of Act in our 1.4 mm con- 
tinuum maps, corresponding to ~ 4mJybearrT 1 , M ~ 1 
and N(H2) ~ 6 x 10 23 crrT 2 in the Southern (proto)cluster, and 
~ 3.2 mJy beam" 1 , M ~ 0.8 Ma 1 and N(H 2 ) ~ 4.5 x 10 23 crrT 2 
in the Northern (proto)cluster. In the latter we detect 7 sources, 
while in the former there are 19 sources detected. Figure[T]shows 
the 1 .4 mm and 3 mm continuum maps for both (proto)clusters, 
with the 26 sources detected at 1.4 mm marked with triangles 
and the two sources only detected at 3 mm marked with squares. 

The properties of the sources are summarized in Table 
Columns 2 and 3 give their absolute positions, the measured 
peak flux intensity and integrated flux density are given in 
columns 4 and 5 respectively. For the unresolved sources, their 
angular size is tipically less than the synthesized beam at the 
respective wavelength. 

Because the brightness temperature at 1.4 mm of the corre- 
sponding Planck function for the strongest source in the region 
is about 2 K, just ~ 2% of the typical hot core temperatures of 
~100 K, we can assume that the emission comes from optically 
thin dust and thus calculate the masses and c olumn densities 
with the approach out l ined b y lHildebrandl (fl983) and adapted by 
iBeuther et all (l2002bl 120051) . We adopted a distance of 2.2 kpc, 
and used a grain emissivity index (3 = 2, corresponding to a dust 
opacity per unit mass K\,\ mm ~ 0.3 and „„„ ~ 0.08 cm 2 g -1 for a 
median grain size a — 0. 1 fim, a grain mass density p — 3 g cm -3 , 
and a gas-to-dust ratio of 186 dDraine et al 1 1200% . 

The selected value of P = 2 for the grain emiss ivity in- 
dex is the value assu med by IBeuther & Schilkel ( 12004 see also 
IBeuther et al.l 120051). and is the typical value found for the 
ISM dust (e.g. lHildebrandlll983HOssenkoDf&Henningll 19941: 
iRodmann et alj|2006l) . Generally. p is thought to range from 1 to 
2. When an IR-to-mm SED is available, as is the case now with 
Herschel results, authors have found p values i n that range and 
even h igher. Examples are /3 ~ 1 .2 and 1 .5 in Ikathborne et al.l 
( 12008b , or p 



1 .7 - 3 3 in lAnderson et al.l d2010b and p 



1.0-1.8 in iRodon et all (l2010h . For IRAS 19410+2336. there 
is no data available to directly calculate the dust opacity index. 
Therefore, in this work we adopted p — 2. 

In contrast to the bolometer observations of the single dish 
data (see below), the interferometric continuum data are pro- 
duced from the line-free part of the spectrum, therefore line 
contribution to the emission is negligible. Potential contribution 
from free-free emission is insignificant as well, since the flux at 
cm wavelengths is on the 1 mJy level and only for the central 
source i n the south, with no emission detected from the other 
sources ( IBeuther et al.l 120031) . Therefore, we can safely assume 
that the mm continuum from which the masses are calculated is 
not severely contaminated by free-free emission. 

The calculated masses and IT column densities are in 
columns 7 and 8 of Table [2] respectively. For their calculation 
we assumed the temperatures shown in column 1 1 of Table [3] 
Those temperatures were determined from the measured H2CO 
line ratios, following the procedure explained and discussed in 
Section l4~2l The visual extinctions in column 9 were cal culated 
assuming A v = N(H 2 )/0.94 x 10 21 dFrerking et al.Lll982l) . 



The masses obtained are strongly affected by the spatial fil- 
tering inherent to interferometers. While varying dust properties 
may account for some of the discrepancy between the masses 
calculated at 1 .4 mm and 3 mm, the spatial filtering is most likely 
the main cause. Comparison with single-dish continuum obser- 
vations of IRAS 19410+2336 at 1.2 mm with the IRAM 30m 
Telescope dBeuther et all l2002bl) show that in our 1 .4 mm con- 
tinuum map of the southern (proto)cluster we are recovering 
only ~ 6% of the flux estimated with single dish observations. 
On the other hand, at 3 mm we recover a relatively larger amount 
of the single dish flux, ~ 37%. This shows that there is a large 
fraction of the mass contained in extended structures, that we are 
filtering out. 

The different percentage of flux filtered out in each wave- 
length band can be explained in the same way as the large (up 
to a factor ~ 10) difference between the calculated masses at 
1 .4 mm and 3 mm. Both wavelengths were observed at the same 
time with the same interferometer configurations meaning that 
although the ground baselines are the same, the uv coverage 
measured in units of wavelength (kA) at 3 mm is more compact 
than at 1.4 mm, therefore tracing more extended components in 
the region. 

In the following, we assume that the missing flux is smoothly 
distributed and affects all cores in the same way. This assumption 
comprises one important caveat when deriving a CMF at high- 
spatial resolution, described in more detail in section \5A\ 



2.18226 10 
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2.182 24 1 5 2.18222 1 5 



E 0.5 




For T k ~ 40 K (see Secs.l3~T1and l4~r2l for more details). 



Velocity (km/s) 



Fig. 2. Comparison of the interferometer (PdBI, filled) and 
single-dish (IRAM 30m, solid line) H 2 CO (3 ,3 - 2 , 2 ) line 
emission towards IRAS 19410+2336. Clearly seen is the "self- 
absorption"-like feature in the PdBI emission, denoting the miss- 
ing flux coming from more extended spatial-scales and filtered 
out by the interferometer. The single-dish spectrum has been 
scaled down for comparison by a factor 10. 



3.2. Formaldehyde 

With the purpose of estimating the temperature structure of 
IRAS 19410+2336, we observed three H 2 CO transitions known 
to be usable as a gas thermometer (e.g., iMangum & Woottenl 
[19931 see Sec. 14.2b . 

Since the PdBI data are clearly affected by missing short 
spacings (an example of this is shown in Figure 0, the region 
was also observed with the IRAM 30m telescope. Figure[3]shows 
the integrated intensity maps of the combined PdBI+IRAM 30m 
data. 
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Fig. 1. Continuum maps of IRAS 19410+2336 obtained with the PdBI. In the top row are the 1 .4 mm maps of the Northern (panel a) 
and Southern (panel c) (proto)clusters. Similarly, in the bottom row are the 3 mm maps of the Northern and Southern (proto)clusters 
in panels b and d respectively. The contouring starts at the 4<x level in all the panels. Contours in panel a increase in lcr steps up to 
8cr and in 4<x steps afterwards, in panel b in 4cr steps, in panel c in lcr steps up to 1 lcr and in 4cr steps afterwards, and in panel d in 
lcr steps up to 16cr and in 4cr steps afterwards. See Table|2]for the cr values. The triangles mark the position of the sources detected 
at 1 .4 mm while the squares are the sources detected at 3 mm. A square appearing in a 1 .4 mm map indicates a source that is only 
detected at 3 mm. Similarly, a triangle in a 3 mm map signals a source detected at the same position in a 1.4 mm map. 



Of the three detected H 2 CO lines, H 2 CO (3 ,3 - 2 o) and 
H2CO (32,2 _ 22,1 ) are the ones with the best signal-to-noise ra- 
tio. According to iMangum & Woottenl (1 19931) those two transi- 
tions are enough to determine temperatures, therefore we did not 
use the third detected transition, H2CO (32,i - 22,0) , for the tem- 
perature analysis. 



Both H 2 CO (3o, 3 - 2 0>2 ) and H 2 CO(32,2 - 2 2 ,i) have their 
strongest emission peak towards the brightest mm source 
detected in the continuum in both the north and south 
(proto)clusters. The secondary peaks in all the maps are in the 
same spatial region as the continuum emission, however there is 
no spatial correlation between the secondary peaks in the H2CO 
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Table 2. Properties of the continuum mm sources in IRAS 19410+2336. 



Source 


R.A. 
(J2000) 


Dec. 
(J2000) 


I v S v 
(mJybeanT 1 ) (mJy) 


T k " 
(K) 


Mass 
(Mo) 


N(H 2 ) 
(10 23 cm" 2 ) 


A v 
(10 3 mag) 


Flags 


1.4 mm with a 1.2" x 0.8" beam 


1-s ... 


19 43 11.553 


23 44 14.15 


4.8 


35 


1.5 


8.2 


0.9 


u 


2-s . . . 


19 43 11.093 


23 44 09.90 


5.7 


35 


1.7 


9.6 


1.0 


u 


3-s ... 


19 43 10.989 


23 44 08.64 


8.1 8.3 


35 


2.5 


14 


1.5 




4-s . . . 


19 43 11.387 


23 44 07.29 


5.3 


35 


1.6 


9.0 


1.0 


u 


5-s ... 


19 43 11.286 


23 44 07.49 


6.2 


35 


1.9 


10 


1.1 


u 


6-s . . . 


1943 11.182 


23 44 07.78 


1.1 16 


35 


4.9 


19 


2.0 




7-s . . . 


19 43 11.093 


23 44 07.41 


7.6 7.7 


35 


2.4 


13 


1.4 




8-s . . . 


19 43 10.989 


23 44 07.33 


5.6 


35 


1.7 


9.4 


1.0 


u 


9-s . . . 


1943 11.191 


23 44 06.68 


1.0 14 


35 


4.3 


17 


1.8 




10-s . . . 


19 43 11.342 


23 44 05.37 


7.2 11 


35 


3.4 


12 


1.3 




1 1-s ... 


19 43 11.229 


23 44 05.37 


5.6 


35 


1.7 


9.5 


1.0 


II 


12-s ... 


19 43 11.357 


23 44 02.39 


1.3 24 


45 


5.5 


17 


1.8 




13-s ... 


1943 11.191 


23 44 03.00 


4.2 75 


90 


8.1 


25 


2.7 




14-s . . . 


19 43 10.950 


23 44 02.15 


4.6 


35 


1.4 


7.8 


0.8 


11 


15-s ... 


19 43 10.837 


23 44 01.94 


4.8 


35 


1.5 


8.1 


0.9 


11 


16-s ... 


19 43 10.614 


23 44 01.33 


7.3 


35 


2.2 


13 


1.4 


11 


17-s . . . 


19 43 10.385 


23 44 02.15 


5.1 


35 


1.6 


8.7 


1.0 


11 


18-s ... 


19 43 11.539 


23 43 57.94 


4.6 


35 


1.4 


7.8 


0.8 


11 


19-s ... 


19 43 11.458 


23 44 57.17 


5.3 


35 


1.6 


9.0 


1.0 


u 


1-n ... 


19 43 10.801 


23 45 08.26 


3.4 


35 


1.0 


5.8 


0.6 


11 


2-n . . . 


19 43 10.728 


23 45 07.41 


6.1 


35 


1.9 


11 


1.0 


11 


3-n . . . 


19 43 10.743 


23 45 06.51 


4.5 


35 


1.4 


7.6 


0.8 


11 


4-n . . . 


19 43 10.706 


23 44 58.42 


2.7 32 


60 


5.3 


25 


2.7 




5-n . . . 


19 43 10.636 


23 44 59.48 


3.6 


50 


0.7 


4.1 


0.4 


11 


6-n . . . 


19 43 10.524 


23 44 58.45 


4.7 


35 


1.4 


7.9 


0.8 


11 


7-n . . . 


19 43 10.469 


23 45 00.2 


3.6 


35 


1.1 


6.1 


0.6 


11 


3 mm with a 2.2" x 1.6" beam 


1-s ... 


19 43 11.558 


23 44 14.41 


5.1 


40 


19.9 


32 


3.4 


11 


20-s . . . 


19 43 11.073 


23 44 12.23 


3.1 


40 


12.1 


19 


2.1 


u, n 


c3mm-sl . . . 


19 43 10.986 


23 44 08.67 


4.8 6.7 


40 


26.2 


29 


3.1 


c(2,3,8-s) 


c3mm-sll . . . 


19 43 11.170 


23 44 07.64 


6.6 13 


40 


50.8 


40 


4.2 


c(5,6,7,9-s) 


c3mm-sIII . . . 


19 43 11.206 


23 44 03.29 


17.7 48 


40 


187.5 


102 


10.8 


c(10,ll,12,13-s) 


16-s ... 


19 43 10.587 


23 44 01.71 


2.9 


40 


11.2 


18 


1.9 


u 


17-s ... 


19 43 10.382 


23 44 02.20 


3.5 3.9 


40 


15.2 


22 


2.3 




21-s ... 


19 43 12.066 


23 44 05.76 


2.0 


40 


7.9 


12 


1.3 


u, n 


22-s . . . 


19 43 12.037 


23 44 03.68 


2.0 


40 


7.7 


12 


1.3 


u, n 


c3mm-nl . . . 


19 43 10.746 


23 45 07.63 


3.2 


40 


12.7 


20 


2.1 


u, c(l,2,3-n) 


8-n . . . 


19 43 10.746 


23 45 05.01 


1.9 


40 


7.5 


12 


1.2 


u, n 


c3mm-nll . . . 


19 43 10.707 


23 44 58.54 


5.1 7.4 


40 


29.0 


31 


3.3 


c(4,5-n) 


6-n . . . 


19 43 10.573 


23 44 58.47 


2.6 


40 


10.2 


16 


1.7 


u 


7-n . . . 


19 43 10.480 


23 45 00.61 


2.2 


40 


8.5 


13 


1.4 


u 



Notes. The last column indicates unresolved cores (u) and detections at 3 mm without a 1.4 mm counterpart (n). The c(X,Y,Z) flag indicates that 
the source is the unresolved combination of sources X, Y and Z. If the source is unresolved, its mass in Column 7 is calculated with the peak 
intensity value of Column 4. 

ia) Derived from the ITCO emission, the 3 mm temperatures are the average in the region (see Sec. 14.1.21 . The values for the 1.4 mm cores are the 
same as in TablefJ] and are repeated here just for reading convenience. 



and dust emission. This prevents us to obtain a temperature for 
each individual core. 

There is an emission feature seen in the northern cluster in 
both H 2 CO (3 ,3 - 2 , 2 ) and H 2 CO (3 2 ,2 - 2 2j i) transitions ~ 4" 
southeast of the main peak, marked with a star in the upper pan- 
els of Figure [3] The feature is at the rest velocity, and we do not 



detect continuum emission neither at 1 .4 mm nor at 3 mm at that 
position. 

The first-moment maps of the combined PdBI+30m data 
(see Fig.|H show only a small velocity dispersion of ~ 1 kms -1 
in both (proto)clusters for the observed H 2 CO lines. Similarly, 
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RA offset [arcsec] 



Fig. 3. Combined PdBI and IRAM 30m integrated emission of H 2 CO towards IRAS 19410+2336. In the left panels are the 
H2CO (3o,3 - 20,2) transition in the Northern (upper panel) and Southern (lower panel) (proto)clusters, and in the right panels 
of the image are the H2CO (32,2 - 22,i) transition. The contour levels are in 10% steps of the peak integrated intensity for each map. 
The triangles mark the 1 .4 mm sources detected, and the 1 .6" x 1 " beam appears in the lower-right corner of each panel. 



the second-moment maps show no strong variation over the 
(proto)clusters. 



3.3. Methyl Cyanide 

We observed the K = - 3 components of the CH 3 CN(6tf - 5^) 
ladder, which are a useful thermometer f or the dense gas (e.g., 
lLoren & Mundvlll984tlZnang et al.l[l998h . We detected CH 3 CN 
emission only towards sources 13-s, 6/7/9-s and 4-n. Figure [5] 
shows the CH3CN spectra toward those three positions. 



The K-0 and K-l components were detected at the three 
positions, and towards 13-s we also detect the K-2 and K=3 
components. Since that is the only position at which those tran- 
sitions are seen, this implies that 13-s is the warmest in the re- 
gion. Also there is a tentative detection of the K-2 component 
towards 4-n. 



Unlike with the H2CO transitions, CH3CN is a dense gas 
tracer, and it is barely affected by the missing short-spacings. 
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Fig. 4. Combined PdBI and IRAM 30m first moment maps of H 2 CO (3 , 3 - 2 0)2 ) towards IRAS 19410+2336 for the Northern (left 
panel) and Southern (right panel) (proto)clusters. The triangles mark the continuum sources detected, and the beam appears in the 
lower-right corner of each panel. There is no signature of a strong velocity dispersion. 



4. Deriving the Core Mass Function 

4.1. Temperature determination 



One of the major caveats of iBeuther & Schilkd (I2004I) in deriv- 
ing a CMF for IRAS 19410+2336 was the assumption of a uni- 
form temperature for both north and south (proto)clusters when 
deriving the masses of the cores. In this work, we try to derive 
a temperature structure for the two subclusters to obtain more 
accurate core masses and by that a more reliable CMF. 

In this section we describe the process used to determine the 
temperature of the detected cores. However it has to be taken 
into account that we are assuming that the kinetic temperatures 
derived from combined PdBI+30m molecular-line analysis cor- 
respond to the dust temperatures of PdBI continuum peaks. As 
is further discussed in Sec. I5.4l this is not ideal, but is one of the 
limitations that we still have to face when working with high- 
spatial resolution mm observations. While combining single- 
dish and interferometer line data is relatively easy, the vastly dif- 
ferent continuum bandpasses of interferometers and single-dish 
bolometers make the combination of continuum data extremely 
difficult. 

4.1.1. Applying H 2 CO line ratios 

Formaldehyde (HhCO), one of the first polyatomic 
molecules discovered in space, has proven its useful- 
ness to derive physical p ropert i es of the interstel la r gas 
(e.g., iMangum & Woottenl Il993t iJansen et all 1 1994 Il995t 



iMiihle et all l2007t IWatanabe & Mitchelll 120081) . In particular, 
H2CO allows one to estimate the kinetic temperature and the 
spatial density in star-forming regions. 

Because H2CO is almost a symmetric rotor (k = -0.96), 
transitions between energy levels with different K are only col- 
lisionally excited. Therefore the comparison between the level 
populations from different K components from the same AJ — 1 
transitio n gives a measure of the k inetic temperature of the 
medium (Ma n gum & WoottenLll993l) . 



In this 
H 2 CO(3 ,3 
ergies E 



up 



work, we used the line intensity ratios of the 

- 20,2) and H2CO (32,2 — 22.1) transitions, with en- 

- 21 K and E 



up 



68 K, respectively. These tran- 
sitions meet the criteria previously mentioned and are only 
~ 254 MHz apart, therefore they could be observed within the 
same spectral setup of the PdBI receivers and the HERA hetero- 
dyne receiver. In this way, possible instrumental uncertainties 
such as telescope efficiency as a function of wavelength, or re- 
ceiver and pointing instabilities, can be disregarded. Hence, one 
can obtain line ratios that are not affected by instrumental effects. 

Since there is only a poor spatial correlation between the 
combined H2CO and continuum emission, we cannot extract 
a temperature for each individual core but only a temperature 
structure of the (proto)clusters as a whole. Only for the brightest 
continuum cores it was possible to associate the H2CO and dust 
emission. Their spectra are shown in Figure |6] 

Having this temperature information is already an improve- 
ment than just assuming a single temperature for all cores. 



7 



J. A. Rodon et al.: Fragmentation in the Massive Star-Fonning Region IRAS 19410+2336 
Table 3. H2CO line parameters for the detected cores. 







H 2 CO(3 


3 _ ^0.2 ; 








2,2 _ z 2, 


) 








Source 


T b 


J T b dv 


v c 


Width 


T b 




V c 




WJ', Ail* 

Width 


Ratio 


T k c 




(K) 


(Kkm s J 


(km s ) 


(kms -1 ) 


(K) 


(Kkms ) 


(km s 


) 


(km s~ 1 ) 




(K) 


1-s" 










0.4 ± 0.2 


1 ± 0.4 


20.4 ± 


0.6 


2.7 ± 1.2 




35 


2-s 


14.7 ± 0.5 


45.7 ± 0.8 


22.2 ±0.1 


3.0 ± 0.1 


1.7 ± 0.3 


4.5 ± 0.8 


21.7 ± 


0.2 


2.5 ± 0.6 


8.5 


35 


3-s 


13.6 ± 0.5 


40.8 ± 0.9 


22.3 ±0.1 


2.8 ±0.1 


1.9 ± 0.4 


3.5 ± 0.4 


21.3 ± 


0. 1 


1.7 ± 0.3 


7.0 


35 


4-s 


17.1 ± 0.8 


51.8 ± 0.8 


22.3 ±0.1 


2.8 ±0.1 


3.0 ± 0.3 


7.4 ± 0.7 


21.9 ± 


0. 1 


2.3 ± 0.3 


5.8 


35 


5-s 


17.5 ± 0.7 


54.4 ± 0.8 


22.3 ±0.1 


2.9 ± 0.1 


2.7 ± 0.2 


7.4 ± 0.8 


21.5 ± 


0. 1 


2.6 ± 0.4 


6.5 


35 


6-s 


18.0 ± 0.9 


56.3 ± 1.1 


22.2 ±0.1 


2.9 ± 0.1 


1.8 ± 0.4 


6.8 ± 0.8 


21.4 ± 


0.2 


3.5 ± 0.6 


9.8 


35 


7-s 


18.1 ± 0.7 


54.6 ± 0.8 


22.2 ±0.1 


2.8 ±0.1 


2.3 ± 0.2 


5.9 ± 0.8 


21.0 ± 


0.2 


2.4 ± 0.4 


7.8 


35 


8-s 


12.4 ± 0.4 


39.0 ± 0.9 


22.3 ±0.1 


3.0 ± 0.1 


1.8 ± 0.3 


4.1 ± 0.4 


20.9 ± 


0.1 


2.2 ± 0.2 


7.0 


35 


9-s 


20.5 + 1.1 


62.7 ± 0.6 


22.2 ±0.1 


2.9 ± 0.1 


2.6 ± 0.5 


7.3 ± 0.5 


21.6 ± 


0. 1 


2.6 ± 0.2 


7.7 


35 


10-s 


18.4 ± 1.3 


61.2 ± 1.3 


22.3 ±0.1 


3.1 ±0.1 


2.7 ± 0.4 


9.3 ± 1.2 


21.5 ± 


0.2 


3.2 ± 0.5 


6.8 


35 


11-s 


19.2 ± 0.7 


58.3 ± 1.0 


22.1 ± 0.1 


2.8 ± 0.1 


3.0 ± 0.5 


8.4 ± 0.9 


21.3 ± 


0.1 


2.6 ± 0.3 


6.5 


35 


12-s 


12.8 ± 0.6 


38.7 ± 0.6 


22.1 ±0.1 


2.8 ± 0.1 


2.6 ± 0.4 


8.0 ± 1.1 


21.3 ± 


0.2 


2.9 ± 0.5 


4.9 


45 


Ij-S 


14. y ± l.U 


oi.j ± U.o 


o 1 a _i_ a 1 
zl.o ± U.l 


1 O _i_ a 1 

J.o ± U. 1 


4.0 ± U.J 


o 1 1 _i_ a q 
Zl. / ± U.o 


ZU.O ± 


u. 1 


i.U ± U. 1 


J. 1 


OA 

yu 


14-s 


5.8 ± 0.2 


13.3 ± 0.7 


22.1 ±0.1 


2.1 ±0.1 


1.7 ± 0.5 


5.1 ± 0.8 


21.2 ± 


0.2 


2.9 ± 0.6 


3.5 


35 


15-s 


6.5 ± 0.7 


21.6 ± 1.1 


22.3 ±0.1 


3.1 ± 0.2 


1.8 ± 0.1 




21.4 ± 


0.1 


unresolved 


3.6 


35 


16-s 


3.5 ± 0.2 


8.6 ± 0.8 


22.2 ± 0.1 


2.3 ± 0.3 


0.8 ± 0.1 




21.3 ± 


0.2 


unresolved 


4.2 


35 


1 IS 


a 1 1 n 1 
4.1 ± U.3 


1 a 1 a n 
10.8 ± U. / 


TO A 1 a 1 

22. 4 ± 0. 1 


2. J ± U.2 


1 1 n 1 
1.0 ± U.l 




22.2 ± 


0.1 


unresolved 


2.3 


33 


lo-S 


a 1 a n 
5.4 ± U.l 


20.3 ± U.9 


A 1 A 1 

22.4 ± 0. 1 


2.0 ± U. 1 


3.1 ± U.l 




21.5 ± 


0.1 


unresolved 


2. / 


35 


19-s 


8.0 ± 0.6 


21.6 ± 0.8 


22.4 ±0.1 


2.5 ±0.1 


1.0 ± 0.3 


3.1 ± 0.5 


22.0 ± 


0.3 


2.8 ± 0.5 


7.9 


35 


1-n 


0.8 ± 0.3 


3.3 ± 0.7 


21.3 ± 0.5 


4.0 ± 1.1 


0.2 ± 0.1 


2.3 ± 0.9 


23.9 ± 1.8 


9.3 ± 3.5 


3.4 


35 


2-n 


1.4 ± 0.2 


4.2 + 0.6 


21.1 ±0.2 


3.0 ± 0.6 


0.8 ±0.1 


1.3 ±0.6 


20.7 ± 


0.3 


1.6 ± 1.1 


1.7 


35 


3-n 


1.4 ± 0.1 


5.1 ±0.7 


21.3 ±0.2 


3.4 ± 0.6 


0.5 ± 0.3 


2.5 ± 0.9 


22.8 ± 


0.9 


4.9 ± 2.1 


2.9 


35 


4-n 


7.8 ± 0.7 


35.6+ 1.1 


21.0 ±0.1 


4.3 ± 0.2 


2.7 ± 0.4 


12.2 ± 1.2 


20.5 ± 


0.2 


4.2 ± 0.5 


2.9 


60 


5-n 


5.4 ± 0.2 


13.6 + 0.9 


21.2 ±0.1 


2.4 ± 0.2 


1.8 ±0.3 


6.0 ± 1.0 


20.2 ± 


0.2 


3.2 ± 0.7 


3.0 


50 


6-n 


6.8 ± 0.2 


20.2 ± 0.9 


21.3 ±0.1 


2.8 ± 0.1 


1.4 ±0.1 


3.4 ± 0.8 


20.7 ± 


0.3 


2.3 ± 0.7 


4.9 


35 


7-n 


5.6 ± 0.7 


16.4 ± 1.3 


20.6 ±0.1 


2.7 ± 0.3 


1.4 ±0.3 


5.9 ± 0.7 


20.5 ± 


0.2 


4.1 ± 0.6 


4.2 


35 



Notes . ( "> No detection in FbCO (3 03 



2o,2) (i) Obtained with a 3" beam to avoid the optically thick regime. <c) From the LVG analysis in section 



4.1.2. LVG modeling 

We extracted spectra towards the positions of each of the 26 
sources detected at 1.4 mm and processed them with CLASS90, 
the obtained line parameters are listed in Table [3] Once we ob- 
tained the R = H 2 CO (3 ,3 - 2 , 2 )/H 2 CO (3 2 , 2 - 2 2 ,i) line ratios, 
they were compared with LVG model predictions of the behav- 
ior of R as a function of molecular hydrogen density n (H 2 ), 
formaldehyde column density N(H 2 CO), and kinetic tempera- 
ture Tk. 

The H 2 density could be estimated from our data. With the 
continuum emission measured towards the positions of the dif- 
ferent sources, we estimate a beam-averaged H 2 column density 
N(H 2 )ranging from ~ 5 x 10 23 crrT 2 to ~ 6 x 10 24 cirT 2 with a 
mean value of ~ 10 24 crrT 2 . For this calc ulation we assumed 
a prio ri a uniform temperature T^ = 46 K dBeufher & Schilkel 
12004 . 

Considering spherical symmetry for the (proto)clusters, at 
the given distance of 2.2 kpc and if is the angular diameter of 
the (proto)cluster, the relationship 



n(H 2 ) = 1.5 x 10 - 16 cirT 



N(H 2 ) 



cm 



kpc 



(1) 



gives an average value n (H 2 ) ~ 10 7 cm 3 for each (proto)cluster. 
For the calculation we adopted an approximate size of 6 ~ 9" for 



the protoclusters, considering the extension of the main 1.4 mm 
emission. This translates to a size of ~ 20000 AU at 2.2 kpc. 

In contrast to n (H 2 ), we cannot directly measure the value of 
N(H 2 CO) from our data. Values for the H 2 CO abundance rela- 
tive to H 2 are reported in the r ange of X (H ? CO) ~ 10~ 9 - 10~ 12 , 
for n ( H 2 ) ~ 10 4 ~ 6 cm' 3 (e.g.JWootten et al.ll 19781; iMundv et al.l 
119871: ICarev et al.lll998t Ivan der Tak et al.ll2000l) . therefore cal- 
culations of an average N(H 2 CO) assuming an X (H 2 CO) value 
from the literature would introduce a high degree of uncertainty. 

To constrain the value of N(H 2 CO) we performed a 
least-squares minimization of the H 2 CO (3o,3 - 2o, 2 ) and 
H 2 CO (3 2 2 - 2 2 1) line intensities as a function of n(H 2 ) 
and N(H 2 CO) for 13-s, the strongest source in the region. 
Assuming Tk = 100K (typical for hot cores), we obtain 
N(H 2 CO) ~ 10 14 ' 5±01 cm- 2 and n(H 2 ) ~ 10 6 - 4±10 cm- 3 . as the 
best solution. This density agrees to within lcr with the previ- 
ously calculated n(H 2 ) ~ 10 7 cm -3 , while the column density 
value can be considered as an upper limit since it was derived 
for the brightest source, with the highest H 2 column density. 

We now compare the LVG model predictions of the be- 
havior of the ratio R as a function of N(H 2 CO) and T^ for 
n(H 2 ) ~ 10 7 cirT 3 with the observed R values. Figure [7] shows 
the comparison. The dotted black contours are the modeled R 
values from 1 to 10 by R = 1 steps, while the solid red con- 
tours are the observed R for sources 3-s, 12-s, 13-s and 19-s (for 
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20 40 60 80 

Velocity (km/s) 

Fig. 5. Observed CH3CN spectra towards sources 13-s; 6,7,9- 
s and 4-n. The detected /if-components are marked with ar- 
rows in each case, and the rest velocity with a dotted line. The 
green solid line is the best Gaussian fit for the lines, obtained 
with CLASS, and the resulting line parameters are in Table 2] 
Only for 13-s were the four ^-components detected, indicating 
a warmer environment. 



viewing simplicity, we only plot a few sources here). The ver- 
tical gray dotted line mark the previously obtained upper limit 
value log N(H 2 CO) ~ 14.5 

For n (H 2 ) ~ 10 7 cm -3 , the calculated upper limit N(H 2 CO) 
puts the cores close to the optically thin/thick regime turnover. 
This turnover is hinted at in the behavior of the contour lines in 
Figure|7] and is located at N(H 2 CO) ~ 10 15 ctrr 2 . At higher val- 
ues the temperature is no longer sensitive to the column density, 
which indicates the onset of the optically thick regime. Another 
method to (qualitatively) estimate the optical depth is by com- 
paring the kinetic temperature Tk with the brightness temper- 
ature Tb. We find that Tb is systematically lower than Tk (see 
Table [3]), therefore if we assume beam-filling the medium is op- 
tically thin. 




Velocity (km/s) 

Fig. 6. H 2 CO spectra of the main sources in each 
(proto)cluster of IRAS 19410+2336. In the upper row are 
the H 2 CO (3 2 2 - 2 2 j) spectra for sources 4-n (left panel) and 
13-s (right panel), and similarly in the bottom row are the 
H 2 CO (3o,3 - 2o j2 ) spectra for those two sources. A Gaussian 
profile was fit to each spectrum, and is shown by a green solid 
line. The dotted line marks the rest velocity adopted for the 
observations. 



Near to the turnover, the uncertainty in the temperature be- 
comes larger, and in this case we see that it is approximately 
~ 15 K. To derive this value we considered the uncertainty of 
cr ~ 0.2 dex in the derivation of N(H 2 CO). Shifting our up- 
per value for N(H 2 CO) by that amount, the derived temperature 
varies by ~ 15 K for 13-s, the brightest source. This can be seen 
in Figure|7] where it can also be noticed that this effect is smaller 
for the sources with lower column density. 

As noted before, there is only a poor correlation between 
the continuum cores and the H 2 CO peaks (except for the main 
sources). Therefore, based on their continuum flux and relative 
location in the (proto)cluster, we separated the cores into sev- 
eral groups. For each of these groups an average value of Tk was 
determined based on the values obtained for each of the sources 
within the group. Column 1 1 of Table[3]contains the kinetic tem- 
perature values assigned to each group which are used to deter- 
mine the masses and N(H 2 ) of the cores shown in Table [2] and 
to build the mass distribution to fit the CMF 

The majority of cores appear to be relatively cold, with a 
temperature of ~ 35 K. Only cores 12-s, 13-s, 4-n and 5-n have 
higher temperatures, up to ~ 90 K in the case of 13-s. This in- 
dicates that these cores are harboring a protostar. In the cases 
of 13-s and 4-n there are reported NIR counterparts (see Sec. 
15.21 ). The average temperature value for the whole region is 
~ 40 + 15 K, in agreement with th e Tk ~ 46 K derived base d 
on IRAS far-infrared observations dBeuther & Schilkel l2004h . 
Consequently, we find that the mean N(H 2 ) does not deviate 
significantly from the value adopted at the beginning despite 
the changes introduced by the different temperature values (see 



9 



J. A. Rodon et al.: Fragmentation in the Massive Star-Fonning Region IRAS 19410+2336 



Table [2J, since that column density was calculated with the 
IRAS-based temperature. 

4.1.3. The CH 3 CN k-ladder spectra 

Methyl cyanide (CH3CN) is a dense gas (n > 10 s cm" 3 ) tracer 
and can also be used as a temperature determinant. As described 
previously in Section 13.31 we only detect CH3CN towards 3 
positions, however only for source 13-s we detect enough k- 
components to derive a temperature, and to that source we apply 
the analysis described in this section. In the other two cases, we 
only detect the first two ^-components, indicating that those are 
cold objects. 

Two methods were used. First we compared the CH3CN 
spectra with LTE m odel spectra produced with XCLASS 
dSchilke et all fl999b . a superset of CLASS of the GILDAS 
package, obtaining a temperature of ~ 80 ± 40 K for 13-s (see 
Fig. [8j. Second, we also de rive a rotational diagram to derive 
the ga s temperature (see e.g jLoren & Mundvll984l;IZhang et"all 
119981) . 

For the rotati onal diagram , we fo llowed the method outlined 
in Appendix B o flZhang et alJdl998l) . Basically, we assume LTE 
and optically thin emission, then the level populations become 
directly proportional to the line intensities of the ^-components 
and are translated into a single temperature via the Boltzmann 
equation. With these assumpt ions, we have the following equa- 
tion (adapted from eq. [B6] of lZhang et al.fl9 98) for the relation 
between the level populations Njjc and the gas temperature T rot 

, Nj, K E lK 1 

In oc — — (2) 

glK k I rot 

where gjjc and k are the statistical weight of the (J, K) level and 
the bo ltzmann constant, respectively. Studies by IWilner et al.l 
(11994 show that T ro i derived with this formulation agrees with 
those obtained with LVG calculations. 

The level population for the upper £-level can be obtained 
from the integrated intensity of the corresponding line. The line 
parameters obtained from the Gaussian fit of the spectra shown 
in Figure|5] and the calculated Nj^k for each level are in Table|4] 
The resulting rotational diagram is shown in Figure|9] along with 
the least-squares fit of all the ^-components using eq. ©, obtain- 
ing the value T rot ~ 100 + 60 K for 13-s. 

From the plot it can be seen that the fit does not represent 
the physical picture well. Explanations of this may be that the 
assumption of optically thin emission is not correct, therefore 
the model applied is not right, or that we are mapping a hot and 
small core with a colder extended envelope, thus presenting two 
different temperature regimes. 

Despite the crude fit, this value agrees very well with the T^ 
obtained from the FLCO LVG model and also with the value 
obtained from the XCLASS model. 

4.2. The Differential Core Mass Function 

Combining the data from both (proto)clusters we derive a differ- 
ential CMF AN /AM, with the number of cores AN per mass bin 
AM. 

One of the strongest caveats (see Sec. 15.4b we faced when 
deriving the CMF was the relatively low number of cores de- 
tected. Because of that, a continuous linear binning in mass was 
not possible. Instead, we performed a logarithmic binning to bet- 



ter represent and analyze the data, meaning that the fixed-width 
mass bins are defined on a logarithmic axis, 

AM = log M k - log M k -i=B (3) 

where B is the constant bin width. Therefore, the k-th mass is 
defined as M k — \0 kB . This binning scheme is further discussed 
in lMafz Apellaniz & Ubedal(l2005l) . 

A priori we do not have a preferred value for B, thus we de- 
rive a mass spectrum for different bin widths, with B ranging 
from 0.001 to 1 in 0.001 steps. Not all the CMFs obtained with 
this method were fitted. To obtain meaningful results, we estab- 
lished the following criteria to be satisfied for a CMF to be fitted: 

- At most only one bin may contain a single core. 

- There must be at least four non-empty bins after the incom- 
pleteness threshold. 

We fit a power-law of the form AN /AM <x M 13 to the CMFs 
satisfying those two criteria, obtaining a Bg index value corre- 
sponding to a given B. The final value B is the weighted mean 
of all the Bb indices with cr < 0.4 and coefficient of correlation^ 
r 2 > 0.9. The different values of Bb satisfying these conditions 

are consistently lower than 2, and have a weigh ted average 

and formal error B = -2.3 + 0.2 (see iRodori 120091 for a more 
detailed description of this method). 

We see then that our results for the CMF slope are steeper 
than the ~ —1.6 found for the Clump Ma ss Functions (e.g., 
IStutzki & Guestedll990tlkramer et al J 19981) . obtained with CO 
emission maps for structures of sizes on the order of ~ 0.1 pc. 
Going to at least one order of magnitude better in spatial res- 
olution, we find that the slope of the Mass Function is consis- 
tently steeper, indicating further fragmentation at smaller spatial 
scales. 

For reference, Figure [l0]shows an example of the CMFs we 
obtained, in this case corresponding to B = 0.214 and having a 
power-law index B = 2.3 + 0.2. A turnover in the distribution 
at the bin centered at ~ 1 .4 M can be seen, containing masses 
above ~ 1 M . It appears in all the derived CMFs at about the 
same position, and since it matches our detection threshold, it 
is likely caused by the low-mass incompleteness of our sample 
rather than being a physical feature. In all of the cases, we fit our 
mass distributions for masses higher than the turnover. 

The incompleteness comes both from both the sensitivity 
limit, as well as the crowding in the region. The sensitivity limit 
has been addressed in Section 13.11 giving a lower complete- 
ness limit of ~ 1 M Q . We find that the effect of crowding can- 
not be properly quantified with the given data. Nevertheless, in 
(sub)mm studies the turnover in the mass distribution and the in- 
co mpleteness of the sample tend to occur at the same position 
(see iKirk et al.ll2006l and references therein), therefore we as- 
sume that our sample is completeness-limited for masses below 
~1M . 

Although we are filtering out considerable fractions of flux 
(section lXTT i on large spatial scales, this filtering likely affects all 
observations in the same way. Although we cannot draw conclu- 
sions about the absolute masses of the cores it is expected that 
their relative masses are correct. 



4 This coefficient is defined as r 2 = 1 - H 22 - , where 5 5 err is the sum 
of squared errors or residual sum of squares, and SS W , is the total sum 
of squares. By this definition r 2 r anges between and 1. A value of 1 
means a perfect fit to the data (e.g, Draper 1998). 
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Table 4. CH3CN line parameters and LTE column densities. 



Source Line 


J T b dv a 

^JV Kill S ) 


iV cr 

^KIIl S ) 


Width cr 

V K1I1 S ) 


T b cr 

tAJ 


logA^ a 
v cm ) 


k = 

13-s 

k = 2 
k = 3 


6.4 1.0 
2.9 0.9 
2.7 0.4 
4.4 0.6 


21.3 0.4 
20.0 0.3 
20.2 0.3 
19.8 0.3 


6.1 1.1 
3.4 0.9 

4.2 0.7 
6.0 1.2 


1.0 0.2 
0.8 0.2 
0.6 0.2 
0.7 0.2 


12.1 0.2 
11.8 0.3 
11.8 0.1 
12.1 0.1 


, 7a £=0 

6 ' 7 ' 9 " S k=l 


2.0 0.2 
1.3 0.2 


22.0 0.1 
20.9 0.1 


1.6 0.2 
1.3 0.3 


1.2 0.1 
0.9 0.1 


11.6 0.1 
11.4 0.2 


4-n * = ° 
k= 1 


2.78 0.02 
1.71 0.01 


21.4 0.4 
20.3 0.3 


4.1 0.8 
2.6 0.7 


0.6 0.1 
0.6 0.1 


11.77 0.01 
11.57 0.01 



120 



100 - 




log 10 [N(H 2 C0)] (cm ) 

Fig. 7. Behavior of the R = H2CO(3o,3 - 2o,2)/H2CO (32,2 - 22 ; i) integrated intensity ratio. The dotted contours represent the R 
values from the LVG model for n (H2) ~ 10 7 crrT 3 (see Sec. 14.1.21 and go from 2 to 10 in unit steps. The red solid contours are the 
observed R values, plotted here for cores 3-s, 12-s, 13-s and 19-s for viewing simplicity. The vertical gray dotted line is the N(H2CO) 
upper limit value obtained with the LVG model. The effect of the optically thin/thick turnover is also seen, which introduces the 
largest uncertainty in the Tk estimation. 



4.2.1. The Cumulative Mass Function 

Using the differential CMF in a low-number sample has the dis- 
advantage that it is sensitive to the arbitrariness of binning, as 
we have shown in the previous section and led to the analysis de- 
scribed there. On the other hand, a cumulative CMF is free from 
any problems associated with the binning, and is more suitable 
for low-number samples. 

With our treatment of the CMF described in the previous sec- 
tion, we have taken into account the arbitrariness of the binning. 



However, for comparison and further analysis we also derived 
a cumulative CMF for our sample. Since our differential CMF 
above ~ 1 M Q can be fit by a single power-law of the form 



A7Y R 
AM 

then the cumulative CMF is 
1 



N(> M) oc 



1+/3 



-M 



(4) 



(5) 
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for/3 < — 1. According to lReid & Wilson! (120061) . when using the 
cumulative CMF one should take into account the upper mass 
cutoff M max of the sample, whether it is a real cutoff or a result 
of the finite sampling. In that case, the cumulative CMF takes 
the form 



N(> M) oc 



J_ (m 1+p - M l+ P\ M < M 

\ + p y vl max Lvl J > lvl ^ LVL m 



0, 



M > U„, 



(6) 



for y8 < -1. However, they also state that a "steep power law" 
(j6 = -2.5) will overwhelm the effect that the upper mass cutoff 
introduces in the fitting. From the differential CMF we see that 
we are close to that "steep power law" limit, so we adopt both 
analytical forms given in equations ((5) and © for the cumulative 
CMF. 

Figure QT| shows the cumulative CMF for masses above 
1 M Q . As stated before we fit both the analytical expressions 
given in equations (0 and (0. The results are /3 = -2.4 + 0. 1 and 
P = -2.2 + 0.1, respectively, with the cr values resulting from 
the fitting algorithm. These results are comparable, showing that 
for a steep power law the inclusion of M max in the definition 
of the cumulative CM F does not affect the result significantly 
(lReid& WirsonLl2006h . Since both results are similar, we adopt 
their average value /3 = -2.3 ± 0.2. 

The cumulative CMF shows that the lower-mass objects 
dominate the fit. This is already visible in the differential CMF, 
but is clearer in the cumulative CMF. This is because the lower- 
mass bins are more populated. We see that of the 25 cores with 
masses above 1 M only 6 are more massive than 2.5 M , there- 
fore the lower-mass bins have more statistical weight in the fit- 
ting. 

There is also an inflection in the distribution starting at 
~ 4M . It is not clear if it is a real physical feature of the dis- 
tribution or a product of the uncertainty of its derivation. We 
discuss this in further detail in Section [ 



5. Discussion 

5.1. Comparison with \Beuther & Schilkd \2004) 



iBeuther & Sdiiike1 ( l2004l) detected 24 mm sources at 1 .3 mm, 12 
in each northern and southern cluster, with a detection threshold 
of 3<t ~ 9mJybeairT 1 . Having better sensitivity, we adopted 
a detection threshold of Act ~ 4mJybeatrr 1 at 1.4 mm find- 
ing 26 cores, 19 in the Southern cluster and 7 in the Northern 
(proto)cluster. Despite the slight difference in the total number, 
the general 1 .4 mm high -resolution cluster structure shown in 
IBeuther & Schilkd d2004l) is recovered. 

This highlights the general mapping difficulties with inter- 
ferometers having only a small number of antennas close to the 
detection limit. For example, contrary to expectations, formal 
3<r limits in such maps are not as reliable as one may expect. 
Therefore, here we raised the threshold to 4<x. ALMA, with its 
many antennas, will overcome such problems. 

From the analysis described in Section 14.21 we ob- 
tain a power-law CMF for IRAS 19410+2336 with an in- 
dex /3 = -2.3 + 0.2, in a greement with the previous result of 
IBeuther & Schilkd d2004 . The new result is more reliable than 
the previous one since with the new data presented here for 
IRAS 19410+2336 we were able to estimate a temperature struc- 
ture for the (proto)clusters. Despite the limitations described in 
sections 14.1.21 and 15.41 that allow us to obtain a more accurate 
value for the core masses and thus avoid the previous caveat of 
assigning a uniform temperature to the whole region. 




1-1037 10 1.1038 10 

Rest Frequency (MHz) 



Fig. 8. XCLASS (dashed line) fit to the observed CH 3 CN spec- 
trum towards source 13-s. 



Nevertheless, this similarity between the CMF indices would 
imply that there is not a strong dependence with the tempera- 
ture structure. If we calculate the cumulative CMF of our re- 
gion, but assuming a uniform temperature T = 35 K, the new 
index is f3 — -2.3 + 0.1. In the same way, if T = 40 K then 
/3 = -2.6 ± 0.2. Both values are indistinguishable within the 
errors from our main result. This feature is in part due to tem- 
perature regime. The functional dependence of the mass with the 
temperature implies that if we were to vary the latter from 10 K 
to 20 K, the mass would change by a about a factor 3. But, going 
to higher temperatures, the variation in mass is less pronounced. 
For example, a change from 40 K to 50 K implies a change in 
mass of less than a factor 1.5, decreasing with increasing tem- 
peratures. Our other uncertainties in the temperature and mass 
calculations have more weight than this effect. Nevertheless, this 
suggests that a detailed temperature structure plays a more sig- 
nificant role in the determination of the mass functions of true 
pre-stellar regions, because of their intrinsic lower temperatures. 
The calculation of the mass function of a star-forming region 
before its members become stars is key in determining how the 
parental cloud further fragments, and wether the IMF is set at 
the fragmentation stage or is a later result of further dynamical 
processes (e.g., competitive accretion, mergers, etc.). 

The calculation of the power-law index was done taking into 
account the arbitrariness of the binning at the moment of deriv- 
ing a differential CMF. Also, the cumulative CMF was taken into 
account in the analysis. 

In the case of the cumulative CMF, the result obtained from 
it matches the result from the differential CMF. Although the 
former is better suited than the latter for the analysis of a low- 
number sample, in this case we observe that the fit to the cu- 
mulative CMF represents the lower-mass end of the distribution 
well but not the higher-mass end. However, a slight 10 K in- 
crease in temperature in the higher-mass end, which is within 
the uncertainty in the temperature determination, flattens the in- 
flection seen starting at ~ 4 M , obtaining also a better fitting for 
the distribution while not introducing substantial modifications 
to the fitted parameters. 

5.2. Continuum Sources 

In the Southern (proto)cluster, there are NIR and MIR 
counterparts detected for se veral of our mm sources. 
iMartm-Hernandez et all (120081) and lOiu et al.l (120081) detect 
in the position of 13-s a bright source in the K s filter and 
the 3.6/um,4.5/um, 5.8yura and 8.0/^m Spitzer/IRAC bands. 
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iMartm-Hernandez et all (120081) suggest that the detected NIR 
and MIR e mission is either leaki ng through an outflow-created 
cavity (see iBeuther et al.l 120031) . or that the cavity itself is 
radiating the emission, based on their estimation that the visual 
extinction at that position should be very large and thus should 
not have any detectable NIR or even MIR emission. We confirm 
their estimation, finding that the visual extinction for 13-s is 
A v ~ 2700 (see Table Source 13-s is in an early stage 
of evolution, according to its NIR ex cess and the presenc e 
of H 2 and Class n CH 3 OH masers (IBeuther et all l2002dh . 
however the det e ction of a VLA 3.6 cm source at its position 
dSridharan et all 120021) . suggest the presence of a recently 
ignited protostar that has already formed an Ultracompact or 
Hypercompact Hn region, the detect ed radio emission being 
consistent with an ionizing B2 V star dMartm-Hernandez et all 
l2008tlOiu et al.Ll2008tlPanagiaLll973h. 

In Figure 2 of lMartm-Hernandez et al! I d2008h . MIR emission 
is also seen at the position of the "subcluster" cent ered in sources 
6-s, 7 -s and 9-s (source mm2 in their work and in IBeuther et al.l 
12001 . There is no detected cm counterpart at that position, sug- 
gesting that either none of the sources have ignited or the ionized 
Hn region is still too small so the free-free emission is confined 
and not detectable. Therefore these sources might be at an even 
earlier stage of evolution as 13-s, the equivalent of a Class 
low-mass protostar, heating their circumstellar dust enough to 
be detected at MI R wavelengths. 

Source nr76 o flMartm-Hernandez et al 1 (12008b is located less 
than 1.5" from source 1-s. However, nr76 is detected in the J, 
H and K s bands and not in any of the Spitzer bands, while 1-s 
is marginally detected at 1 .4 mm and is resolved at 3 mm. The 
lack of emission in the Spitzer bands while being well detected 
in the K s band, and the fact that is a relatively strong source at 
3 mm, make us believe that nr76 is not the counterpart of 1-s, but 
a separate source, appearing nearby due to a projection effect. 



24.5 1 1 , , , 1 , , , , , , , , r 
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Fig. 9. Rotational diagram of the observed ^-levels of 
CH3CN(7 = 6 — > 5) towards source 13-s. The solid line is the 
best fit of Eq. [2] corresponding to a rotational temperature of 
100 + 60 K, in agreement with the kinetic temperature obtained 
with the LVG modeling of the H2CO line emission. 



The case is similar also for source nr71. It is located at ~ 1.5" 
from 20-s, and although nr71 shows strong MIR emission, we 
only detect source 20-s at 3 mm and not at 1 .4 mm. Therefore 
we believe that this is also a case of spatial projection. 

This can be expected, since IRAS 19410+2336 is em- 
bedded in a cluster of Yo ung Stellar Objects (YSOs). 
IMartm-Hernandez et all d2008l) detects 116 NIR/MIR sources 
in a ~ 75" X 75 " region centered on IRAS 19410+2336, while 
lOiu et all (120081) detect 46 YSOs ranging from Class to Class 
n protostars within a radius of ~ 2 pc, in addition to the detection 
of over 800 NIR sources in that region. The fact that there are al- 
ready low-mass stars in the outskirts of IRAS 19410+2336 while 
high-mass stars are still forming in it supports the hypothesis 
that t he low-mass stars f orm before their high-mass counterparts 
(e.g jKumar et alj|2006l) . 

Also in the Northern (p roto)cluste r there are identified 
NIR and MIR counterparts. lOiu et alJ d2008h found 5-band 
Spitzer/IRAC and 2MASS H and K s emission towards 4-n, as 
well as Spitzer/IRAC 3.6/i/n, 4.5/zm, 5.8yum and 8.0pm emission 
towards 1/2/3-n. Source 4-n is likely at an early stage of evolu- 
tion, however it is likely the most evolved source in the Northern 
(proto)cluster. It is not detected at cm wavelengths, therefore it 
has either not yet ignited its protostar or the ionized region is still 
too small and the free-free emission is trapped and thus not yet 
detectable. The other NIR/MIR detection in the (proto)cluster is 
towards sources 1/2/3-n, and their NIR/MIR flux is lower than 
that of 4-n, while their respective 1.4 mm masses are similar. 
This would indicate that 4-n is in a more advanced evolutionary 
stage than 1/2/3-n. 

With a luminosity of ~ 10 4 L Q , IRAS 19410+2336 is a high- 
mass star-forming region, but the (proto)stellar content in its 
core is unknown because of high obscuration. Since it shows cm 
emission, most likely it already has a (proto)stellar component, 
ionizing an Ultra- or even a Hypercompact Hn region. While it 
is not possible to directly measure the masses of the (proto) stars 
with mm data (see Sec. 13. lb , it is possible to estimate the masses 




1 Q 1 1 1 1 I 1 1 1 1 1 1 1 1 I I 

1 10 1 
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Fig. 10. Example of the CMFs obtained, corresponding to B = 
0.214 for core masses above 1.4 M . The solid line is the best fit 
of eq. (01i, with a power-law index B = -2.3 + 0.2. 
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M (Solar mass) 

Fig. 11. Cumulative CMF of IRAS 19410+2336 for core masses 
above 1 M . The solid red line and the dashed blue line represent 
the best fit of equations © with a power-law index/? = -2.2+0.1 
and <j3j with /3 = -2.4 ± 0.1, respectively. The green crosses rep- 
resent a slight increase of 10 K in the higher-mass end, and the 
green dash-dotted line is the corresponding new fit of equation 
(0 (see Sec. 15. U . Notice the flattening of the "bump" in the dis- 
tribution, while the slope of the fit practically does not change. 
In all cases we have normalized A^(> M) by the total number of 
cores N to t- 
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Fig. 12. Comparison between the virial masses and gas masses 
of the resolved cores in IRAS 19410+2336. Virial masses are 
not affected by the interferometric flux filtering, which do affect 
the gas masses. The red dots are the actual values obtained, while 
the blue crosses represent the relationship if we take into account 
~ 10% of the flux filtered out by the interferometer. This is a hint 
of the care that must be taken when interpreting interferometric 
flux values. 

are similar and below the (proto)cluster's Jeans length 
Xj ~ 25000 AU, calculated from the equation 



of the circumstellar structure. Our shortest baseline of 20m at 
the given distance of 2.2 kpc corresponds to a spatial scale of 
~ 38000 AU at 1.4 mm, therefore any structure larger than this 
was not detected. Then the low masses that we calculate from the 
1.4 mm data can be attributed to a circumstellar structure, while 
the surrounding envelope likely contributes to the larger masses 
derived from the 3 mm data. 



5.3. Virial and Jeans analysis 

The large uncertainty that the flux filtering introduces in the mass 
calculations can be seen when comparing the gas masses with 
the virial masses. The former are calculated from the interfero- 
metric continuum emission that is affected by the flux filtering, 
while the latter are calculated using the combined interferomet- 
ric and single-dish H2CO data, with the short-spacings correct- 
ing the effect introduced by the flux filtering. The virial and gas 
masses for the resolved cores in IRAS 19410+2336 are com- 
pared in Figure Q~2] At first glance, the cores do not appear to 
be collapsing, however the gas masses derived from the contin- 
uum are lower than the virial masses because of the missing flux 
discussed previously. If we shift the gas masses higher by a fac- 
tor of 2, taking into account only ~ 10% of the missing flux, we 
see that within the uncertainty on the virial masses now all the 
resolved cores are likely collapsing. 

The relative distances between neighboring cores within 
each (proto)cluster range several thousands of AU. These 



Aj =0.19pc 



10K 



1/2 



n(H 2 ) 
10 4 cm- 



-1/2 



(7) 



for an average n(H2) ~ 10 5 cm 3 and an aver age Ty ~ 40 K 
for th e whole large-scale clumps (see Sec. 14.1.21 Stah ler & Pallal 
12001 . Our resolution of ~ 1" at 1 .4 mm corresponds to a spatial 
scale of ~ 2200 AU, which clearly resolves the Jeans length of 
the region. We are thus resolving the fragmentation of the clump- 
scale (proto)clusters far below their Jeans length. This cannot be 
said in the case of the cores, for which our resolution is above 
the Jeans length of the individual cores, in the 700 - 1200AU 
range. 

We are resolving the Jeans length of the clumps by one order 
of magnitude, almost reaching the Jeans length of the individ- 
ual cores, therefore we can safely assume that we are mapping 
the direct progenitors of single stars or multi ple systems at mos t 
(similar to the Trapezium in Orion, see e.g.. lRodon et al.l I l2008h . 
From the H2CO data, we see that there is only a small signature 
of velocity dispersion (see Fig. |4j, which would imply that the 
(proto)clusters are in a stage of weak dynamical evolution. 

5.4. Considerations 

5.4.1. Masses from (sub)mm interferometry 

All interferometer observations underestimate the masses de- 
tected because of the incomplete sampling of the uv-plane. 
However, this filtering is not necessarily constant throughout the 
observed region. For instance, the single-dish observations map 
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the cores as well as their envelope. This gas and dust reservoir 
is not necessarily gravitationally bound to the cores in its total- 
ity, thus only a fraction of it will be accreted by the cores dur- 
ing their evolution. Theory and observations cannot yet estimate 
how much this fraction is, and in addition, we are not sensitive to 
the mass that is already in (proto)stars or in small, optically thick 
cores. Therefore, the underestimation factor we have calculated 
for the masses is an upper-limit value, obtained in the simple 
approximation that all the mass detected with single-dish obser- 
vations is uniformly distributed and contributes to each core in 
the same way. 

Furthermore, the cores might not draw uniformly from the 
envelope during their accretion phase. The competitive accre - 
tion model (e.g.. iBonnell et alj|2004 120071: iPeters et alj|2010h . 
for example, tells that the most massive cores will accrete more 
from the envelope than lower-mass cores, and that it cannot be 
pointed out to which (proto)star(s) some parts of the envelope 
are going to contribute. It also claims that there maybe is the 
problem that the mass reservoir of each (proto)star is not a well 
defined concept. This has then to be added to the problem of the 
spatial filtering from the interferometer. 

All considered, we believe that the core masses we have de- 
termined are the best values that can be obtained with the cur- 
rently available instruments and techniques. New studies with 
ALMA will certainly improve significantly in this field. 

5.4.2. The nature of the CMF 

In an ideal case, the CMF should be calculated with only the 
pre-stellar (starless) cores because if there is already a proto- 
star, the gas mass accreted onto it is not negligible when com- 
pared with its initial core mass. However, in the high-mass case 
the starless phase is very short lived if existing at all (see e.g. 
iMotte et ail 120071) . Also, some of the MSF theories propose 
that non-prestella r cores would not affect the resulting CMF 
dClark et al.ll2008l see Sec. l531for more details). 

We can assume a relationship of the mass function (MF) 
of a given region to its CMF as CMF - f(M) * MF, where 
f(M) is the fraction we are able to observe. If f(M) = const., 
then our data are a fair representation of the CMF. If it increases 
with mass, for example because a larger fraction is resolved out 
for higher-mass cores, then the slope of the MF we measure is 
steeper than the true CMF. If it instead decreases with Mass, be- 
cause a larger fraction is resolved out for low-mass cores, the 
slope measured will be shallower than the true CMF. The shape 
of this f(M) could be estimated by synthetically "observing" 
model CMFs. However, this is out of the scope of this paper. 

Even without flux losses, unresolved substructure and pro- 
jection effects complicate the picture even more. 

Taking into account all these caveats, we consider that using 
a protostellar CMF as the one derived in this work is the best 
approach we can currently follow for a MSF region, though care 
has to be taken when comparing it with the IMF. ALMA with 
its much better uv-coverage and means to measure the complete 
flux will significantly improve our technical capabilities to study 
the CMF in detail. 



5.5. The Core Mass Function of IRAS 19410+2336 

Most of the studies that trace core scales (~ 0.01 pc) are on 
low-mass star forming regions, since the distances involved are 
of a few hundred parsecs and therefore it is possible to re- 
cover such spatial scales with single-dish telescopes. Some ex- 



amples are lMotte etaD (1 19981) and IJohnstone et al.1 (l2000l) . who 
from single-dish sub-mm observations of the p Ophiuchi region 
the authors derive in each case a mass f unction that res e mbles 
the IMF. A similar result is obtained by iKonvves et alJ d2010h 
with Herschel obse rvations (i.e., in th e far-IR spectrum) of the 
Aquila rift (see also I Andre et al.EoiOl) . Going u p to a few solar- 
like m asses ; there are for example the works of IJohnstone et al.l 
(120061) and INutter & Ward-Thompsonl (120071) on Orion. In this 
case, also from single-dish sub-mm observations, they recover a 
mass function similar to the IMF. 

The mass functions of MSF regions have indeed been de- 
rived in many previous studies, however most of them with 
single-dish telescopes that due to the distances involved of 
more than ~ 2kpc, can only recover scales on the order of 
0.1 pc. According to the criteria we follow that correspond to 
clump scales, i.e., structures that may form (proto)clusters and 
not just a single star or a small multiple system. The results 
show mass spectra with indices 13 between ~ -1.2 and ~ -3 
(e.g. iKramer et alJl!998t iKerton et alJl200U iTothill et alj]2002 



[ Shirley et al.ll2003t |Mookeriea et al.ll2oCTReid & Wilsonll2005l 
2006[lBeltran et al.ll2006HMunoz et al.ll2007t llkeda & Kitamural 
2011) 

The case of IRAS 19410+2336 is, to our knowledge, the 
only one where a high-mass star forming region that could be 
mapped down to core scales presented enough fragmentation to 
derive a CMF. For a MSF region, a mass functio n at similar spa- 
tial scales was derived bv lReid & WUsonl(l2005l) forNGC 7538. 
Their smaller condensations have an average size of ~ 0.06 pc, 
in-between clumps and cores. In fact, they derive a mass func- 
tion with a power-law index fi = -2.0 + 0.3, which despite being 
comparable with the stellar IMF, is also halfway between the 
clump mass function value of -1.6 and the CMF value of -2.3 
found in low-mass star forming regions and in this work for a 
MSF region. 

Remarkably, IRAS 19410+2336 is a very rare object for 
showing more than 10 cores when resolved down to a spatial 
scale of several thousand AU. Similar high-spatial resolution 
(sub)mm observations of MSF regions typically resolve only 
a few cores. For example, interferometric PdBI and SMA mm 
mappings of the MSF regions W3IRS 5, IRAS 06058+2138, 
IRAS 06061+2151, IRAS 05345+3157 and AFGL961 resolve 
cores down to spatial scales between ~ 750 and ~ 5000 AU, 
recovering in all cases fewer than 10 cores in each region 
dFontani et al.l 120091: IWilliams et all 120091: iRodon et alJ 120081: 
lRod6nl2009l) . These regions are similar in distance and luminos- 
ity to IRAS 19410+2336, and although the gas mass varies, even 
the regions more massive than IRAS 19410+2336 show much 
less fragmenta tion at similar spatial sc ales. More recently, there 
is the work of iBontemps et al.l d2010l) in the Cygnus X region. 
They resolve 5 clumps down to spatial scales similar to the ones 
we obtained, recovering 23 fragments in total, from which they 
estimate 9 are probable high-mass protostars. Although the total 
number of cores is similar to this work, each clump fragments to 
'only' 7 cores or less. 

It is unclear why IRAS 19410+2336 is highly fragmented. It 
could be an effect of being at a different evolutionary stage, hav- 
ing a different chemical composition or even through perturba- 
tion and/or a different formation processes. Whichever the case, 
this apparent lack of fragmentation is preventing the derivation 
of a CMF for other individual MSF regions, highlighting the rar- 
ity of IRAS 19410+2336. 

In general, observations of low-mass star-forming regions 
suggest that the IMF of lo w-mass stars is determ i ned at early 
evolutionary stages (e.g., iTesti & Sargentl 119981 : IMotte et all 
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119981 12001b lAlves et al J [20071 iKonvves et alJ 120101: s ee vary- 
ing suggestions e.g.. iGoodwin & Kouwenhovenl 2009h . This is 
not so clear in the case of MSF regions, because dynamical pro- 
cesses like competitive accretion and merging of less massive 
protostars can shape the IMF at later evolutiona ry stages (e.g., 
iBonnell et al]|2004 l2007t iBonneU & Batdl2006l) . 

A similarity between the CMF of a massive star forming re- 
gion and the IMF would suggest that the structure and condi- 
tions within the molecular cloud could determine the IMF. This 
is essentially what was found bv lChabrier & Hennebellel(l2010l) 
in a statistical study of the relationship between the CMF and 
the IMF. They find that both functions are correlated, and they 
argue that the IMF would then be defined in the early stages 
of evolution, being only slightly modified by environmental ef- 
fects. The IMF and the CMF would also be similar if the re- 
lationship between the cores and the stars forming from them is 
one-to-one or nearly one-to-one. This kind of relationship is sup- 
ported by theoretical mo dels explaining the shape of the high- 
mass end of the IMF (e.g. JScalo et al.ll998tlPadoan & Nordlundl 
120021) . and by the apparently constant star formation efficiency 
suggested both by theory and observations dMatzner & McKea . 
l20MlAlvesetaill2007T). 

Also, in the IClark et al.l (120081) picture of competitive ac- 
cretion, the shape of the IMF is independent of when compet- 
itive accretion is halted. This would suggest that the IMF is set 
early on in the evolution of the clump, and therefore having non- 
prestellar cores at the moment of deriving a CMF would not af- 
fect the final results. 

Taking into account the caveats discussed in Section l54l our 
derived power-law index /? is similar to the high-mass end of 
the Salpeter IMF, and one would then be tempted to do a di- 
rect comparison between them. However, the CMF we derive is 
not completely pre-stellar, since we cannot confidently extract 
all the (proto)stellar sources. This means that our CMF is "con- 
taminated" with sources that could already be considered when 
calculating the local IMF in the region. Thus we cannot draw 
solid conclusions about the likeness of the CMF of MSF regions 
and the IMF until the former cannot be derived with just pre- 
stellar cores. 

Even then, caution must be taken when doing a direct com- 
parison between the CMF and the IMF. Numerical simulations 
show that although the overall shape of the IMF in the low- 
mass to intermediate-mass regime is robust against different core 
evolution scenarios, further turbulent fragme ntation of the cores 
may c hange the high-mass slope of the IMF dSwift & Williamsl 
I2008h . 



6. Conclusions 

We resolve the two clumps of the MSF region 
IRAS 19410+2336 into 26 cores at 1.4mm, with a spatial 
resolution of ~ 2200 AU. This resolves the Jeans length of the 
clumps, and the relative distances between the cores are similar 
or smaller than the Jeans length corresponding to individual 
cores. Also, the cores show only a marginal signature of velocity 
dispersion, implying that the (proto)clusters are not in a strong 
dynamical evolution. 

The approach of calculating temperatures using H2CO lines 
might not be suitable for studies at high-spatial resolution and 
density like this one. We find that the H2CO and the contin- 
uum emission only correlate for the brightest source in each 
(proto)cluster and not for the fainter sources, therefore pre- 
venting the derivation of the temperature of each single core. 



However, it allows to derive a temperature structure for the 
(proto)clusters as a whole. 

The temperature structure of IRAS 19410+2336 was deter- 
mined from its H2CO and CH3CN emission allowing the esti- 
mation of a CMF. Taking into account the arbitrariness of the 
mass binning when deriving a mass function, we found a CMF 
index /3 = -2.3+0.2 fo r core masses above ~ 1.4 M H , confirming 
the previous results of iBeuther & Schilkel (120041) with increased 
confidence levels. 

It was not possible to determine with enough confidence the 
evolutionary stage of all the sources detected (i.e, whether they 
are pre-, proto-, or stellar sources). Because of this, the CMF we 
derive is not a fully pre-stellar one, and therefore it cannot be 
considered as a true predecessor of the stellar IMF, but a step 
in between a pre-stellar CMF and the IMF. An additional impor- 
tant caveat is the difficulty to combine continuum single-dish and 
interferometer data because of the vastly different bandpasses. 
This induces missing flux problems severely affecting the mass 
estimates. 

To the present day, mapping of the CMF has been done 
mostly in low-mass star forming regions. We have shown and 
discussed the caveats that involve such a study in MSF regions, 
which in turn explains the rarity of such studies so far. Mapping 
a MSF region down to ~ 0.01 pc scales requires an angular reso- 
lution on the order of one arc-second for the closest regions, and 
tenths of arc-second if we go to larger distances. At (sub)mm 
wavelengths that can only be achieved with interferometry. The 
interferometers operating at this moment normally achieve ~ 1" 
resolution, but to go below it is not that common. Sensitivity 
can become also an issue, and interferometers naturally filter 
out a large part of the incoming flux. When fully operational, 
ALMA will improve significantly on these observational lim- 
itations. Despite the difficulties and assumptions that have to 
be made due to these limitations, the study of the Core Mass 
Function of high-mass star forming regions has great usefulness 
for the study of the origin of the IMF. It is already seen that 
for low-mass stars, the IMF and the CMF are indistinguishable, 
suggesting the kind of relationship between them. That is yet un- 
known for high-mass stars and therefore it needs to be addressed, 
both theoretically and observationally. 
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